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Abstract 



Eshelby showed that if an inclusion is of elliptic or ellipsoidal shape then for any 
uniform elastic loading the field inside the inclusion is uniform. He then conjectured 
that the converse is true, i.e. that if the field inside an inclusion is uniform for all 
uniform loadings, then the inclusion is of elliptic or ellipsoidal shape. We call this the 
weak Eshelby conjecture. In this paper we prove this conjecture in three dimensions. 
In two dimensions, a stronger conjecture, which we call the strong Eshelby conjecture, 
has been proved: If the field inside an inclusion is uniform for a single uniform load- 
ing, then the inclusion is of elliptic shape. We give an alternative proof of Eshelby's 
conjecture in two dimensions using a hodographic transformation. As a consequence 
of the weak Eshelby's conjecture, we prove in two and three dimensions a conjecture 
of Polya and Szego on the isoperimetric inequalities for the polarization tensors. The 
Polya-Szego conjecture asserts that the inclusion whose electrical polarization tensor 
has the minimal trace takes the shape of a disk or a ball. 

Keywords: Polarization tensor, Isoperimetric inequality, Polya and Szego conjecture, Eshelby' conjecture, 
Layer potential 

1 Introduction 

It is well known that amongst all inclusions occupying a given volume the sphere is the 
unique inclusion with maximum surface area. This raises the question as to whether the 
sphere is uniquely optimal with respect to other properties, such as electrical properties. 
It was conjectured by Polya and Szego |48| that the sphere would be the unique inclusion 
minimizing the trace of the electrical polarization tensor when the inclusion and matrix have 
isotropic electrical properties. Here we prove this conjecture. 

A closely related conjecture is the Eshelby conjecture. It is connected to two prob- 
lems. In the transformation problem a region (the "inclusion") in a homogeneous medium 
undergoes a temperature change or phase change which in the absence of the confining 
surrounding medium (the "matrix") would lead to a uniform strain Eq. In the elastic po- 
larization problem, the inclusion has different moduli to that of the matrix, and a uniform 
stress is applied at infinity. If one of these problems has been solved and the field in the 
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inclusion is uniform then (by subtraction or addition of a constant field) one immediately 
has a solution to the other problem. Eshelby |14U15| . following earlier work in special cases 
by Mindlin and Cooper and Robinson [5U], showed that the stress in the inclusion was 
uniform for ellipsoids and furthermore stated (without proof) that "among closed surfaces 
the ellipsoid alone has this convenient property" . Here, for an isotropic matrix, we prove 
Eshelby's conjecture (more precisely what we call the "weak Eshelby conjecture") that the 
inclusion is necessarily ellipsoidal if the field in the inclusion is uniform for all transformation 
strains in the transformation problem, or equivalently if the field in the inclusion is uniform 
for all uniform loadings. In the transformation problem the strain e(x) inside the inclusion 
depends linearly on £o so that we may write e(x) = S(x)eo where S(x) is the fourth order 
Eshelby tensor field. The weak Eshelby conjecture states that if S(x) is constant inside the 
inclusion then the inclusion is an ellipsoid. It follows from certain "trace properties" of the 
second derivative of the Green's function associated with the problem (see, for example, 
equation (6.30) in 36 ) that the isotropic part of S(x) is always uniform and independent 
of the shape of the inclusion |61| . Consequently for any inclusion with a sufficiently high 
degree of symmetry the value of S(x) at its center and the average of S(x) over the inclusion 
equal the value of S in a sphere (or circle in two-dimensions) ^] 123 E] ■ 

For planar elasticity Sendeckyj [HEj, and for antiplane elasticity Ru and Schiavone |52| . 
proved a stronger conjecture (what we call the "strong Eshelby conjecture") that the inclu- 
sion is necessarily elliptical (ellipsoidal) if the field in the inclusion is uniform for a single 
transformation strain in the transformation problem, or equivalently for a single uniform 
loading in the elastic polarization problem. 

Eshelby's conjecture drew increased attention when it was claimed (see 0U1 EU ana - 
references therein) that the field was uniform inside star-shaped polygonal inclusions in 
contradiction to the proof of Sendeckyj. Rodin [51] proved directly that the field cannot 
be uniform inside polygons or polyhedra, and exact expressions for these non-uniform fields 
were later obtained ^] 1271 . Markenscoff showed that the field cannot be uniform if 
any portion of the boundary was planar and that the only small perturbations of any 
ellipsoid boundary that preserve field uniformity in the interior are those which perturb 
the ellipsoid into another ellipsoid Lubarda and Markenscoff 30; showed that the 

field cannot be uniform for inclusions bounded by polynomial surfaces of higher than second 
degree, nor for inclusions bounded by segments of two or more different surfaces, and argued 
that non-convex inclusions are also excluded. 

We remark, in passing, that not only is the field in an ellipsoid uniform for uniform 
loadings but it is also polynomial for polynomial loadings. This was proved for an isotropic 
matrix by Eshelby |15j and for an anisotropic matrix independently by Willis in an unpub- 
lished essay jHU], and by Asaro and Barnett 0. 

Let us now put these conjectures in a precise mathematical framework. Consider in 
R d , d = 2,3 an inclusion f2, which is a bounded Lipschitz domain being inserted into 
a homogeneous medium of conductivity 1 in which there existed a uniform electric field 
E = —a. We assume that the conductivity of fl is k ^ 1. The insertion of the inclusion 
perturbs the uniform electric field and the perturbed electric field is given by E = —Vu 
where the potential u is the solution to 



where a is a constant vector in R d indicating the direction of the uniform field and x(^) 
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Figure 1: The Hashin-Shtrikman bounds for the PTs in R 2 . 



denotes the indicator function of f2. The solution u to Ijl.lf) has a multipole asymptotic 
expansion at infinity, with the leading term being the dipolar one: 



Here u>d is the area of the d — 1 dimensional unit sphere and M is a constant d x d matrix 
independent of a and x. The matrix M = M(fi) := (Afy) is called the polarization tensor 
associated with the inclusion fi. See |51l3*7]. 

In their book 0Hj Polya and Szego conjectured that the inclusion whose polarization 
tensor (PT) has the minimal trace take the shape of a disk or a ball. The purpose of this 
paper is to prove this conjecture in two and three dimensions. In fact, we prove a theorem 
much stronger than the Polya and Szego conjecture. 

In connection with the Polya and Szego conjecture various kinds of isoperimetric inequal- 
ities for the PT have been obtained. See, for example, |46l 1471 155] . The optimal isoperimetric 
inequalities for the PT have been obtained by Lipton j2H], and later by Capdeboscq-Vogelius 
[Hj based on the variational argument in f2B]. The bounds are called the Hashin-Shtrikman 
bounds after names of the scientists who first found the optimal bounds on the effective 
conductivity of isotropic two-phase composites [21], since as pointed out in the caption of 
Figure 2 of [35] , the PT bounds for isotropic M can be obtained as the low volume fraction 
limit of their bounds; more generally for non-isotropic M the PT bounds can be obtained as 
the low volume fraction limit of the bounds of Lurie and Cherkaev E2] and Murat and 
Tartar 021- The PT bounds are given as follows: Let |fi| denote the volume of tt. Then 



u{x) = a ■ x + 



1 (a,Mx) 
~d \x\ d 



+ 0(\x 



), as \x\ — > oo. 



(1.2) 



Tr(M) < |n|(fc- 1+ -), 



(1.3) 



and 



IOITi^M" 1 ) < 



d-l + k 
k - 1 



(1.4) 



where Tr denotes the trace. 

In this paper we prove the following theorem. 



Theorem 1.1 Let be a simply connected bounded Lipschitz domain in R d , d — 2,3. If 
the polarization tensor M(Q) of f2 satisfies the equality in jl.4\) , then Q must be an ellipse 
or an ellipsoid. 

Observe that if a PT M has a minimal trace, then M attains the equality in (|1.4fl and 

M= I (1.5) 

where / is the d x d identity matrix assuming that the volume |f2| = 1. In fact, it can be 
seen clearly from Figure^ which is taken from 3 . In that figure, the horizontal and vertical 
axis represent the eigenvalues of the PT in two dimensions, and hence the constant trace 
lines are those with slope —1. Thus the minimal trace occurs at the unique tangent point 
of the lower hyperbola and a line with slope — 1. This point is an eigenvalue pair of the 
PT associated with the disk. The same argument works for three dimensional case as well. 
Therefore, as an immediate consequence of Theorem 11.11 we obtain the following corollary. 

Corollary 1.2 (The Polya and Szego conjecture) Letil be a simply connected bounded 
Lipschitz domain in M. d , d = 2,3. If 

TrM(Q) = min TrM(D), (1.6) 

where M{D) is the polarization tensor for the domain D and the minimum is taken over all 
the domain with Lipschitz boundary (simply connected or not) with the same volume as £1, 
then is a disk or a ball. 



The concept of the polarization tensor appears in various contexts such as the theory of 
composites (see |37| and references therein) and the study of potential flow 02]. Another 
important usage of the concept is for the inverse boundary value problem to detect diametri- 
cally small inclusions by means of boundary measurements. In fact, one can approximately 
detect, by boundary measurements, the location and the polarization tensor of the inclusion. 
Since the polarization tensor carries important geometric information, such as the volume of 
the inclusion, we are able to recover that information from boundary measurements. It was 
Friedman and Vogelius j^j who first used the polarization tensor for the detection of small 
inclusions. We refer to [5] and references therein for recent developments of this theory. It 
is worthwhile mentioning that the method works for detection of multiple closely spaced 
inclusions 

The main step in proving Theorem II. II is the following theorem from |24j . 

Theorem 1.3 Let be a bounded Lipschitz domain in M. d , d = 2,3. If the polarization 
tensor M(Cl) of fl satisfies the equality in \l-4\) , then for any vector a € M. d the solution u 
to M.l\) is linear in Q. 

Thanks to Theorem 1 1.31 Theorem II .11 is now an immediate consequence of the following 
theorem. 

Theorem 1.4 Let be a simply connected bounded Lipschitz domain in Wi d , d = 2,3. The 
solution u to is linear in Q for any vector a if and only if fl is an ellipse or ellipsoid. 
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It should be noted that only the three dimensional case in Theorem II .41 is new. In two 
dimensions Ru and Schiavone |52| proved a stronger theorem using conformal mappings: If 
the gradient of the solution to l|l.lfl is constant for a single non-zero direction a, then SI is an 
ellipse. This is the anti-plane elasticity case of the strong Eshelby conjecture for elasticity, 
which we now explain in the context of two and three dimensional elasticity. Consider an 
elastic inclusion SI, whose Lame parameters are A, Jl, embedded in a medium in Wi d with 
Lame parameters A, fi. In ^1], Eshelby showed that if SI is an ellipse or an ellipsoid, then 
for any given uniform loading the elastic field inside Q is uniform, and in |15| he conjectured 
that ellipses and ellipsoids are the only domains with this property, which is called Eshelby 's 
uniformity property. 

In order to explain Eshelby's conjecture more precisely, let C = (Cijki) be the elasticity 
tensor of the inclusion-matrix composite, namely, 

C m :=(Ax(K d \H) + Ax(n))tfij<*fci + (MX(R d \fi) +£x(Q))(<^ + 8 a 5 ]k ) . (1.7) 

It is always assumed that 

H>0, rfA + 2/i>0, J1>0 and dX + 2J1 > , (1.8) 

and for technical reasons we also assume that 

(A-A)(/i-#>0, (1.9) 

which means A — A and /i — Jl have the same signs. For given constants d x d matrix A, 
consider the following problem for the Lame system of the linear elasticity: 

V- (C(Vu + Vu T ))) =0 inM d , 

V _ (1.10) 

u(x) — Ax = OQx] 1 d ) as \x\ — > oo. 

If u is the solution to (|1.10(l . then Vu represents the field perturbed due to the presence of 
the inclusion SI under the uniform loading given by V(-Ax). The conductivity model Ijl.lll in 
two dimensions can be regarded as the anti-plane elasticity model of l|1.10f> . What we call 
the strong Eshelby conjecture asserts that if the solution u to l|1.10fl for a single nonzero A is 
linear inside SI, then SI is an ellipse or an ellipsoid. What we call the weak Eshelby conjecture 
states that if the solution u to i|1.10|) is linear inside SI for all A, then SI is an ellipse or an 
ellipsoid. Theorem 1 1 . 41 can be regarded as a solution to the weak Eshelby conjecture for the 
conductivity model. We prove the weak Eshelby conjecture for elasticity. We only state the 
theorem in three dimensions: 

Theorem 1.5 (Weak Eshelby's conjecture in 3D) Let SI be a simply connected bounded 
Lipschitz domain in M 3 . The solution u to \1.1U\I is linear in SI for all A if and only if fl is 
an ellipsoid. 

The strong Eshelby conjecture in two dimensions was proved by Sendeckyj for elasticity 
. In this paper we give a proof of the Eshelby conjecture in two dimensions which is 
completely different from that in |5S]. The novelty of our proof is the use of the hodographic 
transformation. The same approach enables us to construct multiple inclusions satisfying 
Eshelby's uniformity property |23j. 
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Theorem 1.6 (Strong Eshelby conjecture in 2D) Suppose that d = 2. Let il be a 

simply connected bounded domain with the Lipschtz boundary. If the solution u to M.1U\) is 
linear inside fl for a single nonzero A, then Q must be an ellipse. 

We also give an alternative proof of Ru and Schiavone's theorem for the conductivity 
model: 

Theorem 1.7 Suppose that d — 2. Let £1 be a simply connected bounded Lipschitz domain. 
If the solution u to M.l\) for a single vector a ^ is linear in fl, then fl must be an ellipse. 

Our proof uses a hodographic transformation. These have been widely used to solve 
free boundary problems in various problems in mechanics and fluid dynamics, to name one, 
the Saffman- Taylor fingering problem |541 1531 It is also appropriate to mention the 
Vidergauz microstructure. Vigdergauz considered a periodic array of inclusions occupying 
a given volume fraction and found the inclusion shape with minimal overall elastic energy 
|581 1ST?] under certain loadings. He used the fact that the shape would be optimal for one 
of these loadings if the field inside the inclusion was uniform and hydrostatic. See J5| for a 
somewhat simpler treatment. Grabovsky and Kohn proved in the latter paper that the low 
volume fraction limit of the Vigdergauz microstructure is an ellipse. Thus one can expect 
that some variant of Vigdergauz 's complex analytic method might lead us to the proof of the 
strong Eshelby conjecture in two-dimensions. We regard the hodographic transformation 
as such a variant (see also section 23.9 of [37]). Incidentally, we remark that for a dilute 
periodic array of holes under shear loadings the ellipse is not the optimal energy minimizing 
shape [TU] . 

It should be emphasized that the conjectures of Polya-Szego and Eshelby are true only for 
simply connected domains. In a forthcoming paper we construct a family of structures 
with two inclusions in which fields are uniform, and their PT satisfies the lower equality in 

O- 

This paper is organized as follows. In section 2, we review basic facts about single layer 
potentials for harmonic equations and for linear elasticity. In section 3, we prove Theorem 
11.11 We then prove Theorem 1 1.51 in section 3, and Theorem 11.61 and 1 1 . 71 in section 4. 



2 Single layer potentials 

We review some basic facts about single layer potentials for the harmonic equation and 
for isotropic elasticity. We will consider them only in three dimensions. For details of the 
materials presented here, we refer readers to [BJ. 

The single layer potential for the harmonic equation on a bounded Lipschitz domain Q 
in R 3 is defined to be 

Sn[<f>](x) := [ j^da(y), x 6 I 3 , (2.1) 
4tt J dn \x - y\ 

where is a square integrable function on dft and da(y) is the surface measure. Thus Sq[</>] 
is a function on R 3 and <Sfi[0](:r) denotes its value at x. The following boundary behavior 
of the normal derivative of the single layer potential is well-known: 



(,')= (±il + X£)H(.N ;,.(-. ,■ - Oil (2.2) 



± 



G 



where n — {n±, 712, 713) is the outward unit normal to $1, denotes the normal derivative, 
and K.q is defined by 



4tt 



p.v. 



(x - y,n(x)) 
an k — 



<f>(y)dar(y). 



(2.3) 



Here the subscripts + and — denote the limits from the outside and inside f2, respectively, 
and p.v. denotes the Cauchy principal value. See for a proof of i|2.2|) when dfl is smooth 
and [HI] when dfl is Lipschitz. It is known (see also Section 2.47]) that the solution 
u to is given by 

u{x) =a-x + So.[ct>]{x), iel 3 , (2.4) 



where 



Furthermore, we have 



2(^T) 7 ^ 



a • n on 



on 



(2.5) 
(2.6) 



The invertibility of the operator 2 (k\) ^ ~ on L 2 {d£l) is established in |12j . 

It is worthwhile to note that in view of the jump relation l|2.2fl . the usage of the single 
layer potential is natural since l|l.lfl when d = 3 is equivalent to the following problem: 



( Au = 



u\+ 

du 

On 



+ dn 



mnu(R 3 \n), 

on dfl, 
on dSl, 



(2.7) 



, u(x) — a ■ x = 0(\x\ 2 ) 



as \x\ 



00, 



We now review a similar representation formula for isotropic elasticity. The elastostatic 
system corresponding to the Lame constants A, a is defined by 



Cx.^xx := /iAu + (A + /i)V(V • u) . 
The corresponding conormal derivative du/dv on dfl is defined to be 

• \i)n -+- /i(Vu -h Vu )?2 on <9£7 , 



(2.8) 



(2.9) 



where the superscript T denotes the transpose of a matrix. The Kelvin matrix T = (Tij)f - =1 
of the fundamental solution to the Lame system in three dimensions is given by 



r « (z) 



47T tal 



«2 XjZj 
47T Id 3 



x^O, 



where 



1 ( l 

OL\ = - - 



1 



2\[i 2fi + A 



and 0:2 = — ( 



1 



2 V/i 2/1 + A 



(2.10) 
(2.11) 
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The single layer potentials of the density function ip on dfl associated with the Lame pa- 
rameters (A, fi) are defined by 



r - y)<p(y) d <j{y) , x e 



<){> 



Sn[ip](x) 

The single layer potential enjoys the following jump relation: 

= ip on dQ, 



ov 



Ov 



(2.12) 



(2.13) 



where d/dv denotes the conormal derivative defined in (|2.9|l . 

Let \I> be the vector space of all linear solutions of the equation C\^Vl = and du/dv = 
on c?n, or alternatively, 

(2.14) 



f := | ip : dtipj + djipi = 0, 1 < i, j < d 

Here the ipi for i = 1, . . . , d, denote the components of ip. Define 

L%(dfl) := | f € L 2 (dfl) : J f • ip da = for all ip € * j (2.15) 

which is a subspace of codimension 6 in L 2 (dfl). Note that Ijl.lOfl is equivalent to the 
following problem: 



Cx,^u = in W 
C-\ ~u = in fi 



u| _ 
du 



dv 



<9u 

9i> 



on <9il , 
on dfl 

+ 

|-2\ 



(2.16) 



u — Ax ~ 0{\x\ ) as |x| — > oo, 
where £^ ~ and v are the Lame operator and the conormal derivative with respect to the 

Lame constants (A, Ji) of the inclusion. We denote by Sq and Sn the single layer potentials 
on dCl corresponding to the Lame constants (A, /x) and (A, Ji) of the matrix and inclusion, 
respectively. We then have the following representation formula for the solution u to l|1.10f> 
or equivalently 1)2. 16|) : There exists a unique pair (tp, ip) <E L 2 (dCl) x L\ (dCt) such that the 
solution u of H2.16fl is represented by 



u{x) 



' Ax + S n [ip]{x) , xeR 3 \fl, 
^Sn[<p}(x) , xefl, 
where the pair ((p,ip) is the unique solution in L 2 (dfl) x L\{dCl) of 



(2.17) 



SnMl - Sn[ip] 



d_ 

dv 



Si 







--^-Sa[ip] 
ov 



(Ax)\ 9 q 

d{Ax) 



dv 



on on 

on dfl 



(2.18) 



The unique solvability of the integral equation H2.18|l was proved in jlH| under the condition 
(|1.9(l . which is why we assumed this condition. 
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3 Proof of the P61ya-Szeg6 conjecture 



We now prove Theorem ll.4l Theorem l 1.11 follows as an immediate consequence. We consider 
only the three dimensional case because the same proof works for the two dimensional case. 
We begin with the following lemma. 

Lemma 3.1 Suppose that d = 3 and that the solution u to is linear in ft for any 

applied field a G R 3 . Define a linear transformation A on R 3 by A(a) = Vu a |n where u a is 
the solution to Then A is one-to-one and onto. 



Proof. Since the equation in (|l.lfl is linear, that A is linear is obvious. Thus it suffices to 
show that A is one-to-one. Suppose that A(et) = for some a G R 3 . Then Vu Q = in ft. 
Define v(x) = u a (x) — a ■ x for i£l 3 . Then v is the solution to 

f V-(l + (fc-l)x(ft))V<u = inR d , 

\ v(x) = CKM 1 ^) as |a;|->oo, [ " A ] 

By the uniqueness of the solution to (|3.1(l . v = and hence a = 0. This completes the proof. 
□ 

Lemma f3. II can be interpreted as follows: For any vector del 3 there exists a such that 
the solution u to (|1.1|) satisfies 

u(x) = b ■ x + c, ic G ft (3.2) 

for some constant c. It then follows from (|2.4|) and i|2.6|) that 

(b-a)-x + c=(k-l)S n [b-n](x), a: G ft. (3.3) 

In other words, <Sq[6 • n] is linear in ft for any 6 G R 3 . In particular, we have 

S n [nj](x) = linear in ft, j = 1, 2, 3. (3.4) 

By reversing the arguments one can see that (|3.4|) is equivalent to the solution u to (|l.lf> 
being linear in Q for any vector a. 

Now although we have only defined the action of the functional <Sq on scalar functions, 
the obvious generalization of p.lfl defines its action on vector valued functions. In particular 
we have 

5h[n](x) = -V [ — j± :dy, xGft, (3.5) 

Jn - y\ 

which can be seen using the divergence theorem. Thus we get from (|3.4(l 

/ — : -dy — a quadratic polynomial, x G ft. (3-6) 

Jn 47r|a; - y\ 

Since the property (|3.6|) is independent of the conductivity ratio k =/= 1, we have an 
interesting consequence. 

Corollary 3.2 For j — 1,2,3 and < k =/= 1 < oo, let be the solution to f J with 

the 

(k) 

constant in ft for j — 1, 2, 3 and for some k, then Vu) is constant in ft for all k. 



3 

a = ej and conductivity ratio k where {e!,e2,e 3 } is the standard basis for R 3 . // Vuj'' is 
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Suppose that l|3.6|l holds, i.e., there is a symmetric matrix A, a constant vector b, and a 
constant C such that 

1 f 1 

— / -dy = x ■ Ax + b- x + C, x G f2. (3.7) 

4vr 7 n |x - y\ 

This identity has an interesting physical interpretation. If we think of Q as a body of 
constant density in free space then the left hand side is (to within a proportionality factor) 
just the Newtonian gravitational potential at x, and the identity (|3.7I) then says that the 
gradient of this potential, which is the gravitational field, depends linearly on x within the 
body. It has been shown by Dive ^I] an d Nikliborc 0Hj that ellipsoids are the only bodies 
which have this property. Before explaining their proof, let us establish some elementary 
results. 

After a unitary transformation if necessary, we may assume A is diagonal and 
1 f 1 1 3 

— / -dy = - V ajx 2 + b-x + C, xefl, (3.8) 

4vr Jn \x-y\ 2 ^ 3 



for some constants a,j and C, and a constant vector b. We claim that each dj is positive. In 
fact, it follows from 13.511 that 



Sn[rij](x) = —a,jXj — bj, xefl, j = 1,2,3. (3.9) 
Therefore, by l(2~2"j) . we have 

(--/ + K n )[rij] = -djrij on dn, (3.10) 
and hence ^ 

^Wk) 1 -^) [nj]= l + k (k-l)aj nj ° nm - (3 - U) 

Since [—\l + is invertible on Lq(<9£7) := {/ G L 2 (dil) : f dQ f — 0} as was proved in 
|57|. one can see from (|3.10|l that a,j ^ 0. Moreover, the polarization tensor M = (My) 
associated with Q is given by 

Mij = Son Vj {W^T) 1 ~ [ni][y)d(7{v) (3 ' 12) 



M ij = 5 i j\Q,\ - | ^ , i,j = 1,2,3. (3.13) 



as was proved in 0]. Therefore, we have 

k - I 
1 + (fe - l)a 3 ' 

Since the polarization M is positive definite if k > 1 and negative definite if k < 1 regardless 
of fc (see we have cij > 0, j = 1, 2, 3. 

We can now make a complete square out of (|3.8(l and make a translation if necessary to 
conclude that 

3 

13.14) 

i- ' 

with aj > 0. 

The following theorem was proved by Dive ^I] an d Nikliborc 031 for a C 1 domain. The 
same theorem for Lipschitz domains can be proved by a slight variation of their arguments. 



1 f 1 1 ^ 3 , 

— / | -dy = -V ajx] + C, x G 0, 
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Theorem 3.3 Let CI be a bounded domain with a Lipschitz boundary. The relation {3.1$ 
holds if and only if Cl is an ellipsoid of the form 

2 2 2 

fl + | + f§<l. (3.15) 

c l c 2 L 3 

Proof. We briefly sketch the proof. Note that if Cl is an ellipsoid of the form (|3.15|) . then 
(TTH|) holds with 

a, = ™ r ds (3.16) 

for j = 1,2,3 HU IS]. 

To prove the converse, Suppose (|3.14() holds with aj > 0, j = 1,2,3. Then there is a 
unique triple cj., 02,03 satisfying the relation H3.16|l . Existence of such a triple was proved 
in ^202]. Let E be the ellipsoid given by H3.15[) . Then, defining for any region T 



Ny{x) := — / 1 -dy, (3.17) 

4tt J r \x - y\ 



we have 

3 



N E {x) = \^a 3 x] + Ci, xeE, (3.18) 



2 

for some constant C\. For i > 0, let £4 :— {tx\x G -E}. Then by simple scaling one can see 
that 

1 3 _ 

N Et {x) = - a i*j +°t> x£ E u (3.19) 

for some constant Ct depending only on t. Let io be the smallest number such that Cl C E t 
for all t > to. Then there is a point Q which is contained in dE t(l D 9fi, and 

N Et0 \n(x) = N Et0 {x) - JV n (x) = constant, ie!l. (3.20) 

Let n{Q) be the unit outward normal to E ta at Q. Since E to \ Cl lies in one side of the 
tangent plane to E to at Q, we have 

'"^•■W-^Lw*** (3 ' 21 > 

provided that Bj \ 17 is not empty. 

If cWl is C 1 , then n(Q) is also normal to Cl at Q and the normal line goes through Cl. 
Here and in what follows, "goes through Cl" means that there is so > such that the line 
segment {Q - sn(Q)\0 < s < s } C Cl. But by (|3~277jl we get VN E \ n (Q) ■ n(Q) = 0, and 
hence we can conclude that E to \ Cl = by ll3~2"T|) . This is the argument in fTTlllH] . 

If is only Lipschitz, then we can argue as follows. By l|3.21|l . we have VN Et \q(Q) ^ 
provided that E to \Cl is not empty. Thus for any unit vector vo and for any open neighborhood 
V of «o in S 2 , the unit sphere, there is v € V such that 

VN EtgXn (Q)-v^0. (3.22) 
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Choose vq so that the line in the direction vo passing through Q goes through il. Since dfl 
is Lipschitz, there is a neighborhood V in S 2 of vq such that any line in the direction of the 
vector in V passing through Q goes through f2. Then for some v G V i|3.22[l holds. But by 
(|3.2(J|) we get contradiction and hence E to \ O = 0. This completes the proof of Theorem 13. 31 
and hence Theorem 1 1.41 



4 Proof of the weak Eshelby conjecture 

In this section we prove Theorem ll.5l the weak Eshelby conjecture for elasticity. We only will 
prove the three dimensional case because again the same proof works for the two dimensional 
case. 

In fact there is a close link between the weak Eshelby conjecture for elasticity and the 
weak Eshelby conjecture for conductivity (which have just proved). In composite microstruc- 
tures of two isotropic phases it is known that if for some periodic microgeometry the elastic 
field is uniform and hydrostatic (i.e. proportional to the identity) in phase one then that 
microstructure necessarily attains the "bulk modulus type trace bound", equation (6.35) in 
|36| . or the opposite inequality, depending on the moduli of the phases. Then as a conse- 
quence of an argument of Grabovsky [20] (see also section 25.6 of [37|) a solution to the 
conductivity problem (for any direction of the applied field) can be generated from that 
elasticity field, and the electric field in phase one is also necessarily uniform. This argument 
strongly suggests that if Eshelby's uniformity property holds for all applied loadings, then 
there will be one loading for which the field in the inclusion is hydrostatic, and from which 
one can generate solutions to the conductivity problem. As a consequence the electric field 
in the inclusion will be uniform for all applied uniform electric fields, and hence the inclusion 
shape must be ellipsoidal. 

Let us see this directly. Suppose that d ~ 3 and that the solution u to (|1.1Q(I is linear 

in Q for any 3x3 matrix A. Then in the representation formula (|2.17(l . Snip is linear in Q, 

say 

S a [(p](x) =Bx + b, (4.1) 
for some 3x3 matrix B and vector b. It then follows from the integral equation l|2.18|l that 



t av 



= Bx - Ax - 


-b 






on dfl , 


d{Bx) 
, dv 




d{Ax) 
dv 


dU 


on dft . 



(4.2) 



Since C\^Sn[ip] = C\^{Bx — Ax + b) = in f2, the first relation in 14.211 implies that 

Sn{4>}{x) = Bx -Ax + 6, x £ ft. (4.3) 
It then follows from the jump relation (|2.13() 



av 



d(Bx) d(Ba 



dv 



dv 



on. 



Substituting l|4.4() into the first identity in (|4.2() . we arrive at 

'd(Bx) d{Bx) 



di 



dv 



(x) = linear, ieSI. 



(4.4) 



(4.5) 
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The following theorem can be proved in the exactly same manner as Lemma 13. II 



Lemma 4.1 Suppose that d — 3 and that the solution u to $1.101) is linear in Q, for any 
3x3 matrix A. Define a linear transformation A on V , the vector space of all 3x3 matrices, 
by A(A) = Vu^Iq where is the solution to l\l.l(J\) . Then A is one-to-one and onto. 

Since Lemma f4. II means that for any 3x3 matrix B, there is A such that the solution 
u to (|1.10l) satisfies Vu = B in f2, we can conclude that l|4.5|l holds for all 3 x 3 matrix B. 

Let us take a hydrostatic field B, i.e. with B pq 
we have Bx — x. We claim that 



5 pg the identity matrix, in (|4.5|l . Then 



dx 



2^, + 3A 



Sn[nj](x), x £ fi, j = 1,2,3. 



7J 2[i + X 

where the subscript j indicates the j-th component. Similarly we have 

2/z + 3A , 



Sc 



dx 
dv 



2n + A 



■Sn[nj](x), xefl, j = 1,2,3. 



(4.6) 



(4.7) 



We will give proofs of l|4.(j[) and (|4.7|l at the end of this section. 
Because of (14.611 and (|4.7fl . we get 



Sc 



dx dx 
dv dv 



(2^ + 3A)-(2m + 3A) c 
(jc) 3 - = 2^7+A 5 n[niJ(a;), x e 0. 



(4.8) 



Notice that the constant on the right-hand side of (|4.8|l is not zero unless 2/Z + 3A — (2pt + 3A) 
which is excluded by our assumption l|1.9|) . Since the left-hand side of l|4.8|l is linear in ft, 
we deduce that <Sn[nj] is linear in £1 for j = 1, 2, 3. We now conclude from the result in the 
previous section that is an ellipsoid, and the proof is complete. 
Let us now prove l|4.6[l . Since 



— ^ = AV • (x)n + A*(V(x) + V(xf)n = (2/x + 3A)n, 
where n — (ni, n%, n^) is the unit outward normal, we get, for j — 1, 2, 3, 



(4.9) 



dx 
d^ 



(x)j = (2//+3A) 



T(x - y)n{y) 



mi 



3 

(2 M + 3A)V / r je {x-y)n e (y). (4.10) 



Using H2.1()|l we have theorem, 



3 

/ Tji{x -y)n e (y)da{y) 



inn 
ax 



4tt Jan Is - y\ W ^ 



{x - y,n(y)) 
an F — 2/1 



(4.11) 



By Green's theorem, we have 

d 1 



an 



Oj - Vj) 



<9n<, | as - y| 
1 



dn y 



Vj) 



\x - y\ 



\x - y\ 

dy- I [&y(xj - y 3 )\ I z \ dy = 
Ju \ x y\ 



da{y) 



(4.12) 
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(x - y, n(y)) J _ / ,_ A _ f rij{y) 



and hence 

L {x > - ^'^w^' - - 7„ ra^ 1 "'- <413) 

It then follows from (|4.11() and the fact that ot\ — a 2 — 2 ^+\ th & t 

tj B J,^-vW,)My) = Z^J^ My) . -gi^feM (4,4, 

By substituting this back in l|4.10|l we get 14. 6|) . 
To prove l|4.7|l . it suffices to note that 

^ = AV • (x)n + fl(V(x) + V(x) T )n = (2J1 + 3A)ra. (4.15) 
ov 

The proof is now complete. 

5 Alternative Proofs of the strong Eshelby conjecture 
in 2D 

In this section we give alternative proofs of Theorem 11.61 and 11.71 We first prove a lemma 
on univalence of the analytic functions which will be used to prove Theorem 11.61 and 11.71 
The proof of the following lemma relies on the level curve argument which was used in 
various contexts. We particularly mention the work of Alessandrini and Nesi [2] in which 
they showed the univalence of cr-harmonic mappings in the context of periodic composite 
materials. The two phase inclusion-matrix problem in the free space can be viewed as a low 
volume limit of the periodic case. 

Lemma 5.1 Let Q be a simply connected domain with Lipschitz boundary in C and let 
f be a analytic function in C \ fi such that there are constants a =/= and (3 such that 
f(z) — (uz + (3) — > as \z\ — > oo. If f(z) — ix 2 for z — x\ + ix 2 G dfl, then f is univalent 
in C \ fi. 



Proof. Let u and v be the real and imaginary parts of /, respectively. Observe that / maps 
C* \ £1 onto C* \ [ic\,ic2\ where c\ — min z6 gsi v(z) and C2 = max z6 gn v(z), and C* is the 
Riemann sphere. Since C* \ Q and C* \ [ic\ 1 ic2\ are simply connected, it suffices to show 
that f'(z) ^ for z G C* \ fl to prove univalence of /. It is obvious that /'(oo) ^ 0. We 
will show that Vu(z) ^ for any z G C \ f2. 

Suppose that Vm(zo) = for some zo G C \ f2. Then there is an integer m > 2 such that 
u takes the form 

oo 

u{z) — u(zq) = \z - z \ n (a n cos n9 + b n sin nO) , z near zp, (5-1) 

n—m 

where z — zq = \z — zo\e and a 2 m + b 2 m ^ 0. Therefore, there are 2m branches of level curves 
C{zq) :— { z | u(z) = u(zq) } coming out of the point zq. Since u(z) — axi+bx2+c+0(\z\~ 1 ) 
as \z\ — > oo for some real constants a, b, and c, at most two of these 2m branches are 
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unbounded and extend to the infinity. For the other branches of C(zq) which are bounded, 
one of the following two occurs: (i) a branch intersects dfl, (ii) a branch meets another 
branch and makes a loop. 

Suppose that (ii) occurs. Then the loop may encircle or be contained in C \ f2. If it 
is contained in C \ fi, then by the maximum principle, u is constant inside the loop and we 
have contradiction. If the loop encircles f2, then u is constant (zero) on dfl and on the loop 
(the constants may be different). Let O be the annular region enclosed by the loop and d£l. 
Then by the Cauchy-Riemann equations, we have = on 90, and hence v is constant in 
0, which leads us to a contradiction. 

Now suppose that (ii) does not occur for any bounded branch. Let lj, j — l,...,fc 
(k > 2) be bounded branches of C(zq). Let Zj be the point where lj intersects dfl. Then 
there are two branches, say l\ and I2, such that l\, I2, and the arc on d£l connecting z\ and 
z-i make the boundary of a connected region f2. Since u = u(z\) = u(z2) = on l\ and I2, it 
follows from the maximum principle that u = in Q, which contradicts to the assumption 
that a^O. This completes the proof. □ 

We now provide the alternative proofs of Theorem 1 1 . 61 and 1 1 . 71 
Proof of Theorem 11. 71 Let u be the solution to (|1.1J) for some a ^ 0, and suppose that u is 
linear in fl. Put u l := u\q and u e :— u\ R2 \ji- It is proved in [221 that there are functions U l 



and U e analytic in Q and R \ fl, respectively, such that Re U l 



~-U l 



k-1 — 
——U l ■■ 



u e 



on dQ, 



and RcU e = u e , and 



(5.2) 



for some real constant c. Since u l is linear in Q, there are complex numbers 8 and 7 such 
that 

U l {z) = 6z + j. 



Define /n by 



fn(z) := U e (z) 



k+l 



(Sz + j) + ic, zeC\Q. 



Then by (|5.2|) we have 

We now define i/'n by 

i>n{z) 



fn(z) 



- 1 t-l(5z + j), zedfl. 



k — 1 



+ -z, zeC\n. 



S(k - 1) 

Then one can see from (|5.5|) that 

V'fiW — i x 2i for z — xi + 1x2 G dQ. 
Since u e {x) — a ■ x — 0(|a;| _1 ) as \x\ — > 00, ipn takes the form 

ipn(z) = az + + (pa(z), 



(5.3) 
(5.4) 

(5.5) 
(5.6) 
(5.7) 



(5.8) 



for some complex numbers a and (3, and a analytic function <pn in C \ such that <pn(z) 
0{\z\- r ) as \z\ -> 00. 
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We claim that In fact, if a — 0, then h 



= Re ipn is the solution to 



Ah = 
h = 

>(^)-c = 0(|x|- 1 ) 



I 2 
on 951 



iii A' \ O, 



(5.9) 



as \x\ 



for some real constant c. But the maximum principle (with the point at infinity being 
regarded as a point on the Ricmann sphere) implies c = and hence h = in R 2 \ fl. Thus 
we conclude that a =/= 0. 

It then follows from Lemma |5 . II that ipa is a univalent mapping from C* \ onto C* \ 
[ici,ic2] where c\ = min z6 an Irn V'n (z) and ci — max^gn Imi/^z). We may assume that 
ci = — 1 and C2 = 1 by scaling if necessary. 

If D is the unit disk, we can construct such a mapping ipn explicitly. In fact, it is the 
Koebe function 

M*) = \(z-\), 1*1 >1- (5-10) 

Then ifin is a univalent mapping from C* \ D onto C* \ [— i, i]. 
Let 

F(z):=^ 1o ^d(z), \z\>1, (5.11) 

in which i/jq 1 is the hodographic transformation. Then F is a univalent mapping from \z\ > 1 
onto C \ and F(po) — oo. Moreover, F(z) behaves as a linear analytic function at infinity. 
Since C* \ D is simply connected, it follows from Caratheodory's theorem P-18] that F 
extends to C* \ D as a continuous function and F is a homeomorphism from 3D onto dfl. 
Observe that if \z\ = 1, then F(z) = ^^(1x2), and hence ImF(z) = X2- In other words, 
Im F is a solution to 



Au = 

u(x) = X2 

u{x) - (bixi + 6 2 a; 2 + c) = Od^l" 1 ) 



I 2 
on 9Z? 



(5.12) 



as \x\ 



00, 



where 61, 62, and c are real constants. 
We claim that c = and 

ImF(z) = bxxi + b 2 x 2 
In fact, if we define u by 



-61X1 + (1 - 6 2 )x 2 



z > 1. 



w(z) = ImF(z) - 



bixi + b 2 x 2 



-b\Xi + (1 - 62)^2 



\z\ > 1, 



(5.13) 



(5.14) 



then u is a solution to (|5.9|l with replaced by D, and hence c = and it = 0. 
We now get from l|5.13|l that 

F{z)= 1Z + C +-, |*| >1, 



(5.15) 



for some constants 7^0 and [3. One can easily see that the image of the unit disk under F 
is an ellipse. This completes the proof of Theorem 1 1.71 □ 
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Proof of Theorem 11.61 Let u be the solution to (|1.1U|) for some constant ay, not all zero, 
and assume that u is linear in £1. We first invoke the following complex representation of 
the solution to i|1.10[l from (see also [5] Theorem 6.20]): Let u = (it, v) be the solution 
of (|1.1(J|) for d — 2 and let u e :— u| C wj and := u|n- Then there are unique functions tp e 
and ipe analytic in C \ f2 and (fi and tpi analytic in such that 



2fi(u e + iv e )(z) = Kcp e (z) - z(p' e (z) - ip e (z) , z £ C\Q , (5.16) 



2/u,(ui + ivi)(z) = Kifi(z) - zip'^z) - i\)i{z) , zeQ, (5-17) 

where 



A + 3/i _ A + 3/i 
A + ii A + ju 

Moreover, the following holds on dft: 

1 



(5.18) 



— ^K<p e {z) -Z(p' e {z) - ip e (z)J = — ( Ktpi{z) - ztp'^z) - tpi{z) ) , (5.10) 



ip e (z) + ztp' e (z) + ip e (z) = tpi{z) + zip' t (z) + ^i(z) + c , (5.20) 

where c is a constant. Equation (|5.19|l is the continuity of the displacement and l|5.20fl is 
the continuity of the traction. 

It follows from l(5~TO|l and (|5~2U|) that 

(k + l)<p e (z) = (y + + (l - |) (^( z ) + M*)) + c, z£» (5.21) 

Since u is linear in 17, y>j and ^ are linear analytic functions in Q by the uniqueness of tfi 
and ^ij. It then follows from l|5.21fl that there are complex numbers a, (3, C such that 

<p e (z) = az -Pz + C, zedQ. (5.22) 

Suppose that (3^0. Define /n by 

fn(z) = ^[Ve(z)-az-C]+^z, 2 gC\fi. (5.23) 

Then fn(z) — 1x2 for z = xj + 1x2 G <9f2. Following the same argument as in the proof of 
Theorem 1 1.71 we conclude that f2 is an ellipse. 

If j3 = 0, then ip e can be extended to the whole of C as an entire function. In fact, if we 
define (p e (z) = az + C in f2, then tp e is analytic in C \ d£l and continuous on dQ, and hence 
is an entire function (It can be proved using Morera's theorem that that if a function 
defined in an open set U is analytic in U minus a Lipschitz curve and is continuous in U , 
then that function is analytic in U .) But, since u(x) — aijXiej = 0(|x| _1 ) as x — > 00, 
ip e takes the form ip e {z) = 7z + / e (z) for some constant 7 and a analytic function f e in C\f2 
such that / e (z) = 0( 1 2; I —1 ) as |z| — > 00. So / e extends as a bounded entire function and 
hence f e is constant and the constant is 0. Thus 

tp e (z) =~{z for z G C \ O. (5.24) 
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It then follows from (|5.19|) and (|5.20|l that 



7 -z - -^-ipe(z) = ( ntpi(z) - zip\{z) - ipi(z) ) , (5.25) 



2/i 2/i 2/i \ 



(1 + i)z + Mz) = <Pi(z) + z^{z) + Aiz) + c , (5.26) 

on <9f2. Therefore there are complex constant S, n, C such that 

i) e (z) = 5z-rJz + C. (5.27) 

If rj = 0, then by the same reasoning as to derive Q5.24|) . one can see ip'( z ) is constant. 
Therefore, u e + iv e is linear, and hence u(a;) — ■ aijXiej = for all x £ C \ fi. This is 
possible only when f2 is an empty set. Thus, r\ ^ 0. Now let 

gn{z) = ^[^ e {z) - 5z - C] +^z, zeC\Q. (5.28) 

Then gn{z) — 1x2 for z = x± + 1x2 G 917. Following the same argument as in the proof of 
Theorem 1 1.71 we conclude that fl is an ellipse. This completes the proof. □ 

We finally mention that the strong Eshclby conjecture for three dimensions (except in 
the special case where the field in the inclusion is hydrostatic) has not been proven, not even 
for the conductivity case. 
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